Spin waves in zigzag graphene nanoribbons and the stability of edge ferromagnetism 
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We study the low energy spin excitations of zigzag graphene nanoribbons of varying width. We 
find their energy dispersion at small wave vector to be dominated by antiferromagnetic correlations 
between the ribbon's edges, in accrodance with previous calculations. We point out that spin wave 
lifetimes are very long due to the semi-conducting nature of the electrically neutral nanoribbons. 
However, application of very modest gate voltages cause a discontinuous transition to a regime of 
finite spin wave lifetime. By further increasing doping the ferromagnetic alignments along the edge 
become unstable against transverse spin fluctuations. This makes the experimental detection of 
ferromagnetism is this class of systems very delicate, and poses a difficult challenge to the possible 
uses of these nanoribbons as basis for spintronic devices. 



Graphene is being hailed as the big promise for nano- 
electronics and spintronics. Its unique transport proper- 
ties are expected to play a fundamental role in the de- 
velpment of new technologies [l|43| . New physics is also 
emerging from the interplay between low dimensionality, 
a bipartite lattice and electron-electron interaction. One 
of the most striking properties of graphene nanoribbons is 
the possibility of spontaneous magnetization |^-|yl- This, 
combined with the long spin-coherence times of electrons 
propagating across graphene, indicates that this system 
is a strong candidate for future spintronics applications. 

The ground state properties of magnetic graphene 
nanoribbons have been extensively explored by a variety 
of methods. Recent works have investigated the proper- 
ties of static excited states based on adiabatic approx- 
imations [3, Q. This approach has been employed to 
describe the lowest-lying excitations of magnetic metals 
with relative success. However, it is well known that it 
misses important features of the excited states, such as 
its finite lifetime. This arises due to the coupling between 
spin waves and Stoner excitations, a distinctive feature 
of itinerant magnets. Moreover, these recent investiga- 
tions of excited states seem to have disregarded the an- 
tiferromagnetic coupling between the magnetizations on 
opposite edges of graphene nanoribbons. As we shall see, 
this leads to an incorrect prediction concerning the wave 
vector dependence of low energy spin excitations. This 
has already been demonstrated more than a decade ago 
in the seminal work by Wakabayashi et al. |9|. Those 
authors used an itinerant model to describe the tt elec- 
trons in graphene nanoribbons of various widths. They 
showed clearly the presence of a linear term in the spin 
wave dispersion relation for small wave vector. 

One interesting feature of magnetic graphene nanorib- 
bons is that the spins along each border are ferromag- 
netically coupled to each other, but there is an antifer- 
romagnetic exchange coupling between the two opposite 
borders. This coupling is mediated by the conduction 
electrons, and decreases as the ribbon width is increased. 
Thus, it may appear, at first sight, that this antiferro- 
magnetic coupling should be unimportant in wide rib- 



bons. It has been shown, however, that this coupling 
is extremely long ranged in graphene and other related 
materials [10l4l4|. Thus, even in rather wide nanoribbons 



this coupling asserts itself, as we shall see. 

We describe the electrons in graphene using a Hubbard 
model. 
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where tij are hopping integrals (i 7^ j) and on-site en- 
ergies (i = j), U is the effective intra-atomic Coulomb 
interaction and cl^ creates one electron at the atomic 
state at site i with spin a. Here we only consider near- 
est neighbor hoppings. We took [/ = 2 eV (~ 0.77t), 
as in Ref. [15[. This model provides a good qualita- 
tive description of tt electrons in graphene, as well as 
the magnetic effects deriving from the screened Coulomb 
interaction. The magnetic ground state is described self- 
consistently within a mean-field approximation. We im- 
pose local charge neutrality on every atom in the ribbon 
and determine the magnetic moment of each atom in the 
unit cell individually. We find that the magnetic moment 
of the edge atoms is 0.24/xb and decays rapidly towards 
the center of the ribbon, in close agreement with calcula- 
tions based on density functional theory [7|. Notice that 
we consider the effective Coulomb interaction to be ac- 
tive in every atom in the system. The fact that the mag- 
netization is essentially localized at the edges emerges 
naturally from our self-consistent treatment. 

The spin excitations are extracted from the properties 
of the transverse dynamic susceptibility. 
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and S^ = {S^y are the spin rais- 
ing and lowering operators. By treating the Coulomb 
interaction term within a random phase approximation 
we obtain a closed equation of motion for x^ (^) (the 
Fourier transform of X^~Ji}) in terms of the mean-field 
susceptibility x'"^+"(^) [ll, 

x+-{n) = \i + ux^"^+-{n)r\^"^+-{n), (3) 
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FIG. 1: (Color online) Schematic depiction of the zigzag 
nanoribbon's geometry. The dotted lines encircle two arbi- 
trary unit cells, labeled m and m' . The indices I and I' refer 
to atoms inside each unit cell (as in equation (5]) . 



where x^ and x £^rc matrices comprising all mag- 

netic sites in the system and / is the identity matrix with 
the same dimension as x"' and -^ ■ The notation 

used in the equation above is schematic. A more precise 
statement on the form of the calculated susceptibility is 
given below. 

The spectral density 



A,{n) = -3x«(^), 



(4) 



where 9 denotes the imaginary part, may be interpreted 
as the density of states of magnons in the system. The 
dynamic susceptibility just described is the response of 
the system to an externally applied field of frequency Q, 
transverse to the ground state magnetization direction; 
spin waves appear as peaks in the spectral density. 

The graphcne nanoribbons we study have translation 
symmetry along the ribbon length (denoted here by x), 
but not along the ribbon width (y). It is convenient to 
define a mixed Bloch-Wannier basis to describe the elec- 
tronic states, 
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where c{xm,yi) is the annihilation operator for a Wan- 
nier state at a site I in unit cell m, a = VSao is the 
distance between unit cells along the ribbon length and 
Co « 1.42 A is the carbon-carbon distance. In this 
representation, the transverse dynamical susceptibility 
X;"yr(ri; k) is a matrix, where I, V label sites within a unit 
cell; each element of such matrix is a function of the wave 
vector k along the length of the ribbon, as well as of the 
energy i7. The unit cell is depicted in Fig. [1] 

We start by discussing the spectral density Aiik] Vt) for 
a nanoribbon of fixed width. In Fig.[2]we plot the spectral 
density as a function of Vl for fixed values of wave vector 
fc, for a ribbon with 16 atoms in the cross-section. The 
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FIG. 2: (Color online) (a) Spectral density associated with 
spin waves, projected on the "up" edge, for a neutral ribbon 
8 atoms wide, for selected wave vectors (indicated in the fig- 
ure); (b) linewidth as a function of wave vector; (c) spin wave 
dispersion relation deduced from the peaks of the spectral 
density (squares). The dashed curve is a plot of the quadratic 
dispersion relation found in adiabatic calculations 



main contribution for the excitations should come from 
the edges, where most of the magnetic moment of the 
system is concentrated. Thus, we only need to plot the 
spectral density projected on the "up" edge, that we label 
as i = 1 (the spectral density at the "down" edge has 
similar behavior) . Spin wave energies increase with wave 
vector, as usual, but the k dependence is not quadratic, 
as would be expected from a simple fcrromagnet. 

A plot of the dispersion relation deduced from the peak 
positions (Fig. [2][c)) shows that the dispersion is linear 
quite far out into the Brillouin zone (20% of the zone 
boundary), and in fact is quasi linear at large wave vec- 
tors as well. This may be understood if we map the spin 
degrees of freedom of this system onto a simple effective 
spin model, as illustrated in ref. [9|. 

Calculations of spin excitation energies based on an 
adiabatic approach have been reported recently [7| . They 
find a quadratic energy-wave vector dispersion relation 
with a spin wave exchange stiffness of 320 meVA^. We 
plot this dispersion relation in Fig. [2][c) for compari- 
son. Although the energies found using the adiabatic 
approach are of the same order as those obtained via 
dynamical calculations, the discrepancy between the de- 
pendencies on wave vector is remarkable. 

In Fig. [2ja) we see that the spin wave peaks are ex- 
tremely narrow for small wave vectors, indicating a very 
large spin wave lifetime. This is compatible with the ex- 
istence of a threshold for Stoner excitations of the order 
of the semiconducting gap in these ribbons. Only spin 
waves with energies equal to or larger than this gap are 



damped. This means that the low energy spin dynam- 
ics (represented by long wavelength spin waves) may be 
well described by effective localized spins hamiltonians, 
but as the wavelength of the excitation become smaller 
the itinerant nature of the system reveals itself. Thus, 
a simple ferromagnetic Heisenberg hamiltonian is clearly 
not the appropriate model to describe the spin degrees 
of freedom of this fascinating system. 

One very attractive feature of graphenc is the possibil- 
ity of controlling its carrier density by electrostatic gat- 
ing. In the present case this feature opens up a very 
exciting possibility: by controlling the electron density 
we may be able to tune the relaxation time of spin ex- 
citations in graphene. As shown above, long wavelength 
spin waves are essentially undamped in electrically neu- 
tral graphene ribbons. In fig. [3] we show that very modest 
changes in the electron density can induce rather large 
damping, reducing considerably the relaxation time for 
spin excitations and shifting their energies. The ori- 
gin of this damping is simple to grasp: the density of 
Stoner modes is very small at small energies in undoped 
graphene ribbons due to the fact that the density of states 
p near the Fermi level Ep is zero (the antiferromagnetic, 
undoped nanoribbon is semiconducting). As the den- 
sity of states is increased by the gate voltage, p is in- 
creased for energies close to Ep, giving rise to a signif- 
icant enhancement of the density of Stoner modes. As 
it is well known 16l-ll9|. spin wave damping in itiner- 
ant systems occurs through the decaying of magnons into 
Stoner modes, a mechanism very similar to the Landau 
damping of plasmons in metals. In Fig. ^ih) we show 
how the density of Stoner modes at the spin wave ener- 
gies (given by the spectral function A^ associated with 
the non- interacting susceptibility x^'^'^^ ) is enhanced by 
increasing electron density. 

It is also clear from Fig. [3] that the extra damping is 
accompanied by a shift in the spin wave energy. Once 
again, this is related to the enhancement of the den- 
sity of Stoner modes, via the Kramers-Kronig relation. 
The non-interacting susceptibility -^ enters in the 

denominator of the dynamic susceptibility x^ i ^^ ii^" 
dicated in Eq. [3] Its imaginary part is responsible for 
the finite lifetime of spin waves in itinerant magnets; its 
real part produces a shift in the spin wave frequencies, in 
much the same way as dissipative forces shift the natural 
frequency of mechanical oscillators. Thus, enhancement 
of damping also implies a larger frequency shift. 

There is another facet to the onset of strong spin wave 
damping in graphene nanoribbons. Spin excitations with 
infinite (or extremely long) lifetimes are associated with 
strongly localized spins, whereas strongly damped spin 
waves are found in systems where magnetism is itinerant 
in nature. It is very rare that one system can be tuned 
to be either a localized or an itinerant magnet with the 
change of a single parameter, easily accessible experimen- 
tally. It is an extremely exciting prospect that this kind 





5 

°b 1 2 3 4" 
Sn (10 ^ electrons/atom) 



90 100 110 
n (meV) 



20 40 
Q. (meV) 



FIG. 3: (Color online) (a) Spectral densities for q = 0.266 A~^ 
and different doping levels (0, 1, 2, 3 and 4 xlO^'^ elec- 
trons/atom) for a ribbon 8 atoms wide. Larger doping means 
smaller energy and smaller peak height, (b) Linewidth as a 
function of doping (solid circles, scale on the left) and density 
of Stoner modes A^ at the spin wave energies (solid triangles, 
scale on the right), (c) Spectral density at low energies for 
the same wave vector and doping levels as panel (a). Notice 
that the spectral density is absolutely flat in this region for 
zero doping. 



of control is available in graphene nanoribbons. 

The lifetime of spin waves in zigzag graphenc nanorib- 
bons can be dramatically reduced, as we just saw, by very 
modest doping (as small as 10""^ electons/atom). By fur- 
ther increasing doping we notice that the ferromagnetic 
alignment along the borders becomes unstable. A sign of 
this instability is the appearance of a very soft spin wave 
mode as doping increases, as can be seen in Fig. ^ic). 

The instability of the ferromagnetic alignment can be 
confirmed by the behavior of the mean-field transverse 
susceptibility at zero frequency, as a function of wave 
vector, Al{k) = x^°^+"(fc,f^ = 0) Q. In a stable ferro- 
magnetic system AQ(fc) has a single maximum at A: = 0, 
as illustrated by the zero doping curve in figure Sl^a). As 
doping increases, a peak develops close to fc = 0, until, at 
large enough doping (in this case, 0.01 electrons/atom), 
a pronounced maximum appears at a finite value of k. 
The existence of peaks in A^ (fc) at finite values of k means 
that the true ground state of this system is not ferromag- 
netic along the edges, but most probably a spin density 
wave characterized by those finite wave vectors. 

One virtue of our simple model is that we can tune 
parameters and explore various behaviors. By changing 
the strength U of the Coulomb interaction we noticed 
that the doping level at which the instability appears 
changes. This is illustrated in Fig. 2] for two different 
values of U\ ioi U ^ 2 eV the ferromagnetic alignment 
becomes unstable at ~ 5 x 10~^ electrons/atom, whereas 
for [/ = 1.5 eV it is stable for doping levels as large as 
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FIG. 4: (Color online) The zero frequency mean-field spectral 
density Ao(fc) for different doping levels and two values of the 
effective Coulomb interaction, [7 = 2 eV (a) and U = 1.5 eV 
(b). The curves have been displaced vertically for the sake of 
clarity. In panel (a) the doping levels are 0, 2, 3, 4, 5 and 10 
me/atom; in panel (b) the doping levels are 0, 5, 10, 15 and 
20 me/atom. Doping increases from bottom to top. 



1.5 X 10"^ electrons/atom. It would be interesting to 
build a U X dn phase diagram, but our intention here is 
to point out the dependence and the general trend. 

The stability analysis we performed is complementary 
to that presented in references [15|, |2l|, [22 1 , where energy 
differences between coUinear and non-collinear configura- 
tions in the direction transverse to the ribbon were anal- 
ized. To the best of our knowledge, ours is the first anal- 
ysis that take into account the possibility of non-collinear 
ordering along the ribbon's edges. 

We have investigated wider nanoribbons (up to 32 
atoms wide, although this is only limited by computer 
time) . The most important effects of increasing the rib- 
bon width are i) the enhancement of a quadratic con- 
tribution to the spin wave dispersion relation (due to 
the partial suppression of antiferromagnetic coupling be- 
tween edges) and ii) the reduction of the doping range 
(for fixed U) within which ferromagnetism along the bor- 
ders is stable. 

In conclusion, we have shown that the lifetimes of spin 
excitations in zigzag graphene nanoribbons can be tuned 
by the application of modest gate voltages. This allows, 
at least in principle, electrical control of the magnetic re- 
laxation rate. We have also demonstrated that there is a 
sharp transition between the character of the spin excita- 
tions in neutral and doped nanoribbons: while in neutral 
ribbons the long wavelength excitations have essentially 
infinite lifetime (a feature shared with magnetic insula- 
tors), any amount of doping, however small, leads to fi- 
nite lifetimes (as in magnetic metals) . Finally, we showed 
that further increasing doping makes the ferromagnetic 
alignment unstable against transverse spin fiuctuations, 
a fact that should be carefully invesigatcd if these sys- 



tems are to be used in technological applications. We 
are confident that our results open very exciting possibil- 
ities both in spintronics technology and for fundamental 
understanding of magnetism at the nanoscale. 
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